Strain tuning of topological band order in cubic semiconductors 
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We theoretically explore the possibility of tuning the topological order of cubic diamond/zinc-blende semi- 
conductors with external strain. Based on the tight-binding model, we analyze the evolution of the cubic semi- 
conductor band structure under hydrostatic or biaxial lattice expansion, by which a generic guiding principle is 
established that lattice expansion can induce a topological phase transition of small band-gap cubic semicon- 
ductors via a band inversion, and further breaking of the cubic symmetry leads to a topological insulating phase. 
Using density functional theory calculations, we demonstrate that a prototype topological trivial semiconductor, 
InSb, is converted to a nontrivial topological semiconductor with a 2% — 3% biaxial lattice expansion. 



PACS numbers: 71.20.Nr, 71.70.Fk, 73.20.At, 77.65.Ly 
I. INTRODUCTION 



The recent discovery of topological insulators has gener- 
ated great interest in the fields of condensed matter physics 
and materials science, largely driven by their exotic surface 
electronic properties. 13 These materials are insulating in the 
bulk, but support topologically protected helical edge or sur- 
face state in two- or three-dimensional systems. 1-5 The sur- 
face states are predicted to have special properties that could 
be useful for practical applications in spintronics and quantum 
computation. 6 To date, topological states have been observed 
in a number of materials, 7-35 opening the door for detailed 
and systematic investigations of various topological phenom- 
ena in laboratories. While the essential features of topologi- 
cal insulators have been established, our understanding of the 
fundamental physics in relation to basic material properties is 
far from complete. This type of knowledge will be crucial in 
potential device applications. 

In this work, we theoretically demonstrate the possibility 
of tuning the topological order in cubic semiconductors by 
applying external strain, using InSb as a prototype system. 
Our choice of the material system is motivated by the fact 
that many technologically important semiconductors crystal- 
ize in one of the cubic crystal structures, including the group- 
IV elements in the diamond structure, and the III-V and II- 
VI compounds in the zinc-blende structure. Ternary com- 
pound semiconductors with the tetrahedral coordination, such 
as half-Heuslers 16 ~ 18 and chalcopyrites 19 , also belong to this 
category. We first establish, based on tight-binding analysis, 
that lattice expansion can induce a topological phase transi- 
tion of small band-gap cubic semiconductors via a band in- 
version, and further breaking of the cubic symmetry leads to a 
topological insulating phase. Using density functional theory 
calculations, we then demonstrate that InSb enters into a topo- 
logically nontrivial state at a reasonable biaxial lattice expan- 
sion around 2% — 3%. The relationship between mechanical 
strain and topological order has been discussed in layered bis- 
muth compounds 36 as well as half-Heuslers 1 . Our analysis 
provide a physical picture of these findings. 



II. GENERAL CONSIDERATIONS 

The band topology of materials with time-reversal symme- 
try is characterized by the so-called Z-i topological index. In 
general, the calculation of the Z 2 index requires the knowl- 
edge of the occupied Bloch wave functions in the entire Bril- 
louin zone because it is a global property of the energy bands. 
However, in most practical situations the index can be de- 
termined by counting the number of band inversions at the 
time-reversal invariant fc -points. In passing from a topologi- 
cally trivial to nontrivial state, a band inversion, i.e., switch- 
ing of occupied and unoccupied bands with opposite parity 
around the Fermi level, must happen. If an odd number of 
band inversions occur within the Brillouin zone, then the ma- 
terial could be a topological insulator. 

We now specialize to zinc -blende semiconductors. In a 
typical band structure, the low-energy electronic properties 
is dominated by bands around the T point. Away from the 
r point the valence and conduction bands are well separated. 
Therefore we only need consider the band order at V. The 
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FIG. 1. (Color online) The band order of cubic diamond/zinc-blende 
semiconductors at the T point with Aso > 0. (a) The normal band 
order, (b)-(c) The inverted band order can be obtained from the nor- 
mal band order by either increasing the spin-orbit coupling Aso (b) 
or by decreasing the coupling potentials V ss and V pp (c). The band 
inversion strength is defined as AE — Er 6 — Er s , which has pos- 
itive value in normal band order and negative value in inverted band 
order. 
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relevant states are the Tq anti-bonding state of the s-orbitals, 
and the r> and Tg bonding states of the p-orbitals. Using the 
tight-binding model, 37 ' 38 the band energies of these states are 
obtained 



(a) static lattice (b) a=a +1%a (c) a=a +2%a (d) a=a +3%a 
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where e s i,2 an d £pi.2 are the s- and p-orbital energies of the 
cation and anion, respectively, V^ s and V pp are the coupling 
potentials of the s-s and p-p bonds, and Aso is the spin-orbit 
coupling strength. Note that and g states have opposite 
parity. 

Usually, the state is located above the Ty^ states, form- 
ing the conduction band minimum. This is the so-called nor- 
mal band order [see Fig. 1(a)]. There are two ways to change 
the band order, (i) If the spin-orbit coupling strength is in- 
creased, depending on the sign of Aso, either the Tg (Aso > 0) 
or the Tj (Aso < 0) state could rise above the Tq state, re- 
alizing the necessary band inversion [see Fig. 1(b)]. (ii) In 
addition to varying Aso, which is an intrinsic property of the 
material, the band order can be also changed by varying the 
coupling potentials V ss and V pp [see Fig. 1(c)]. Increasing 
the lattice constant leads to a decrease of the coupling poten- 
tials, which lowers the T 6 anti-bonding state and raises the 

8 bonding states. At sufficiently large lattice constant, the 
band inversion can be realized. Thus in the design of topolog- 
ical insulators, both factors should be taken into account. 

The relation between the band gap and lattice constant has 
been discussed previously in the context of band gap pressure 
coefficient. 39 ' 40 In a more realistic situation, there are other 
factors that could affect the band positions, such as the p-d 
coupling in the compounds with active d-orbitals in the va- 
lence bands, as well as the kinetic energy effect under pres- 
sure. A detailed analysis can be found in Refs. 39 and 40; the 
general trend we discussed above still holds. 

Finally, if Ep 8 > Ef e in the inverted band structure (Aso > 
0), the twofold degenerate Tg state is fully occupied while the 
fourfold degenerate r§ state is half filled, resulting in a zero 
band-gap semiconductor [for example, see Fig. 2(c) and 2(d)]. 
To realize a truly insulating ground state, a possible way is to 
apply a biaxial strain so as to break the cubic symmetry and 
thereby open a band gap at V point [see Fig. 3(d)]. 



III. METHODOLOGY 

Having established the general tuning principle, next we 
verify whether the anticipated topological phase transition can 
happen in real materials at a reasonable lattice expansion. We 
choose InSb as a prototype system because it has a relatively 
small band gap of 0.235 eV, 48 indicating that InSb is on the 
verge of becoming a topological insulator. 
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FIG. 2. (Color online) The band structures of InSb at the equilibrium 
lattice constant (a) and with 1% (b), 2% (c), and 3% (d) hydrostatic 
lattice expansions. The Te, IV, and Tg state are denoted by red, 
black, and blue color, respectively. The size of dots is proportional 
to the probability of s-orbit projection. 



TABLE I. The relative energy changes of states Er 6 , Er 7 , and Er 8 
under the hydrostatic lattice expansions (a — b — c) ranging from 
0%to4%. 
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The band structure calculations in this work were per- 
formed using full-potential linearized augmented plane-wave 
method 41 ' 42 implemented in package WIEN2K. 43 A converged 
ground state was obtained using 10 4 fc -points in the first Bril- 
louin zone with i?MT^max = 9.0, where Rmt represents the 
smallest muffin-tin radius and K nydx is the maximum size of 
reciprocal-lattice vectors. The muffin-tin radius was set to 2.5 
Bohr for both In and Sb atoms. Wave functions and poten- 
tials inside the atomic sphere are expanded in spherical har- 
monics up to I = 10 and 4, respectively. Spin-orbit cou- 
pling is included by a second-variational procedure, 41 where 
states up to 9 Ry above Fermi level are included in the ba- 
sis expansion. The modified Becke-Johnson exchange po- 
tential together with local-density approximation for the cor- 
relation potential (MBJLDA) 47 was used to obtain the band 
structures. The MBJDA potential can effectively mimic the 
behavior of orbital-dependent potential around the band gap, 
and it is expected to obtain accurate positions of states near 
the band edge 47 . In fact, both standard local-density approx- 
imation (LDA) 44 ' 45 and generalized gradient approximation 
(GGA) 46 calculations predict that InSb is a zero-band gap 
semiconductor with an inverted band order, while the MB- 
JLDA calculation produce a band gap about 0.253 eV, which 
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FIG. 3. (Color online) The evolution of band order of InSb with 
nonhydrostatic expansion (a) in afe-plane and (b) along c-axis. The 
black solid circles are the results without lattice optimization, orange 
down-diamonds and purple up-diamonds are optimized by LDA and 
GGA potential receptively, (c) and (d) are normal band order and 
inverted band order, which correspond to the states marked by dotted 
circles in (a). The size of dots is proportional to the probability of 
s-orbit projection. 



is in excellent agreement with the experimental value of 0.235 
e y 48 rpjjg MBJLAD method has been applied for small band- 
gap half-Heuslers in the context of topological insulators 20 ' 21 
and yields favorable comparison with more sophisticated GW 
method. 49 

We consider two types of lattice expansion. The first is 
the hydrostatic lattice expansion, equally increasing the lat- 
tice constants along all the three axes. This situation is more 
of academic interest rather than for practical applications. In 
the second case, we consider the biaxial lattice expansion, 
where the crystal structure is expanded in the a6-plane and 
relaxed in the c-axis. This can be achieved by growing the 
sample on a substrate with a larger lattice constant. Although 
the MBJLDA potential gives very accurate band gap and band 
order, it has no corresponding definition about its energy func- 
tional. 47 Hence, both the LDA and GGA for the exchange- 
correlation potential were used to optimize the lattice con- 
stants. 



IV. RESULT AND DISCUSSION 

The crystal structure of InSb is zinc -blende with space 
group F43m (No. 216). The unit cell of InSb contains two 
inequivalent atoms, In and Sb, with the fractional coordinates 
(0,0,0) and (0.25,0.25,0.25), respectively. The initial equilib- 
rium lattice constant of 6.479 A is taken from previous exper- 
imental measurements. 48 

Figure 2 shows the band structures of InSb at its equilibrium 
lattice constant and under hydrostatic lattice expansion rang- 
ing from 1% to 3%. The relative energy changes of the states 



Et 6 , Er 7 , and Er 8 are listed in Table I. We find that, as antic- 
ipated, the band inversion takes place upon lattice expansion, 
somewhere between 1% and 2%. In addition, we find that 
the spin-orbit coupling strength, characterized by Er s — Er 7 , 
stays almost a constant with a slight decrease as the lattice 
constant increases. Hence the evolution of the band order dur- 
ing lattice expansion is dominated by the change of the cou- 
pling potentials V^ s and V pp . 

Although a hydrostatic expansion already converts the nor- 
mal band order into the inverted band order, a nonhydrostatic 
strain is still needed to creat a band gap at T point by break- 
ing the cubic symmetry. As shown in Fig. 3, we investigate 
the evolution of band order with the nonhydrostatic expan- 
sion both in afe-plane and along the c-axis. Note that in the 
tetragonal symmetry the band inversion strength AE is rede- 
fined as the energy difference between the s-like states and 
the top of valence bands. First, we consider the lattice expan- 
sion in the a6-plane with fixed lattice constant along the c-axis 
(without lattice optimization). The AE change from positive 
values to negative values as shown in Fig. 3(a) (black solid cir- 
cles). The band structures of two typical states are presented 
in Fig. 3(c) and 3(d). One is normal band order with 1% ex- 
pansion [Fig. 3(c)] and another is inverted band order with 3% 
expansion [Fig. 3(d)]. Then, we consider the lattice optimiza- 
tion along the c-axis by using both LDA and GGA potential. 
It is well known that LDA usually underestimates the lattice 
constant and GGA overestimates it, hence these results give 
us a lower and upper bound. The true transition point for the 
expansion in aft-plane is in the range of 2% ~ 3%. Finally, if 
we expand the lattice along c-axis, there exist similar topolog- 
ical order transition. The only quantitative difference is that a 
larger lattice expansion is needed, as shown in 3(b). 

Finally, we mention the possibility that the tensile biaxial 
expansion may also be realized by growing the cubic semi- 
conductor on the top of a piezoelectric substrate whose lattice 
can be changed by applying an electric field. In this way, the 
topological order of the semiconductor can be tuned with elec- 
tric field. 



V. SUMMARY 

In summary, we have explored the strain tuning of topolog- 
ical band order in cubic semiconductors by the combination 
of simple tight-binding analysis and density functional theory 
calculation. We have predicted that InSb can realize the topo- 
logical insulating state under the biaxial lattice expansion of 
2% — 3%. This work provides a generic guiding principle for 
tuning the topological order of the cubic semiconductors and 
offers an opportunity for experimentally exploring the proper- 
ties of the topological surface states on such technologically 
relevant materials for practical applications. 
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